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ABSTRACT 

We have modelled the light distribution in 10095 galaxies from the Millennium Galaxy- 
Catalogue (MGC), providing publically available structural catalogues for a large, rep- 
resentative sample of galaxies in the local Universe. Three different models were used: 
(1) a single Sersic function for the whole galaxy, (2) a bulge-disc decomposition model 
using a de Vaucouleurs (i? 1 / 4 ) bulge plus exponential disc, (3) a bulge-disc decomposi- 
tion model using a Sersic (i? 1//n ) bulge plus exponential disc. Repeat observations for 
~ 700 galaxies demonstrate that stable measurements can be obtained for object com- 
ponents with a half-light radius comparable to, or larger than, the seeing ZiaZf-width 
at half maximum. We show that with careful quality control, robust measurements 
can be obtained for large samples such as the MGC. We use the catalogues to show 
that the galaxy colour bimodality is due to the two-component nature of galaxies (i.e. 
bulges and discs) and not to two distinct galaxy populations. We conclude that under- 
standing galaxy evolution demands the routine bulge-disc decomposition of the giant 
galaxy population at all redshifts. 

Key words: astronomical data bases: catalogues - galaxies: general - galaxies: fun- 
damental parameters - galaxies: structure - galaxies: statistics. 



1 INTRODUCTION 

To understand the origins of the diverse galaxy population 
observed today it is essential to quantify the properties of 
different structures which may be assoc iated with separate 
formation processes. Recent work (e.g. iBaldrv ct al. 2004) 
reported that the local galaxy population c onsist s of two 
distinct classes (red and blue). iDriver et all J2006T) demon- 
strate that the galaxy bimodality is particularly distinct in 
the colour-log(n) plane, where n is the Sersic index for the 
best-fitting 7? 1//n model to the 2D galaxy light distribution. 
Furthermore, they suggest that the bimodality may be bet- 
ter interpreted as a representation of the two-component na- 
ture of galaxies (i.e. red bulges and blue discs). They arrive 
at this conclusion by noting that E/SO's and Sd/Irr's lie in 
the two distinct peaks, whereas the Sabc's (i.e. bulge+disc 
systems) straddle the two peaks. The colour bimodality may 
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reflect the end-result of two different processes associated 
with bulge and disc formation. 

It has long been known that although galaxies can cover 
an expanse of quite different morphologies, they have many 
prevalent features in common, most notably spheroids (or 
bulges) and discs. The prominence of these features is of- 
ten used to classify gal axies and place them into an evo- 
lutionary scenario ?e.g. lHubblelll925 Ide Vaucouleurslll959l : 
Ivan den Bergbll976T) . Early classification of galaxy morphol- 
ogy was based on the visual inspection of images, and al- 
though this can be useful, significant differe nces between 
individual classifiers hav e been observed (e.g. lLahav et alJ 
ll995tlNaim et alJll995ah . and less subjective techniques are 
clearly desirable. 

In order to provide a more quantitative assessment of 
galaxy morphology, a number of different methods have been 
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its asymmetry in a C — A system |Burbid ge et al 
Elmegreen et alJll992t iDoi et alJ fl993l : ISchade et al 
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phology in a way that is not dependent on an assumed 
model. Extended classification systems also include a mea- 
su re of the 'dum piness' of the galaxy light (C — A — 
S; IConselicdEooli) . Other methods employ artificial neu- 
ral networks (e.g. | Naim et alJ Il995bl : lLahav et al] Il996 



CL996) to classify galaxies according to a 
Hubble sequence scheme. 

Bulge-disc decomposition is another popular and useful 
method for quantifying the morphologies of galaxies by fit- 
ting model surface brightness profiles to data. Galaxies are 
described in terms of the two most easily recognised stellar 
components: bulges and discs. Such an approach has some 
clear disadvantages; the fitting process is model-dependent 
and does not account for secondary features such as spi- 
ral arms, rings, bars, star forming knots, tidal tails, and 
other asymmetries. However, even at z ~ 1, bulges and 
discs are clearly in place and appear to be the dominant 



galaxies (e.g. 


Simard et al.lll999: Ravindranath et al. 


iBarden et al. 


120051 iKoo et al.1 2005 


iTruiillo et al. 



In addition, software is now available which makes 
it fairly straightforward to perform bu lge-disc decomposi- 
tion on many thousands of galaxies (e .g.lSimard et alJl2002t 
[ Peng et al] (20021: ITruiillo fc Aguerril 12004 Ide Souza et al] 

A bulge-disc model is also convenient because the sur- 
face brightness profiles of bulges are known to be wel l mod- 
elled by R 1/4 ide Vaucouleur dll94gft. R 1/n . Sersic JSersid 
|l_968|), or exp onential laws JXiidredakis fc Sanders! Il994 
Ide Jondl996ah and discs are observed to follow a pu re expo- 
nential profile Jde Vaucouleuril959tlFreemanlll97rJl . Bulges 
and discs are also often observed to have different average 
colours, metallicities, and kinematics, justifying their treat- 
ment as distinct entities. Most models of galaxy formation 
and evolution also in volve separate fo rmation scenarios for 
bulges and discs (e.g. lCole et aljEoOOt) . 

Using HST imaging, several statistically significant 
samples of high redshift galaxies with measured structural 
parameters have been constructed (e.g. [Simard et al.l 120021: 
iRavindranath et al jl20ol IBarden et aljboOSft . However, for 
the local Universe there only exist samples of hundreds of 



galaxies, that have had a bulge-disc decomposition, and 
these are often preselected to be only late- type (e.g. Ide Jond 


ll996al:lGraham & de Blokl 


200lHMacArthur et al.ll2003f). or 


early- type (e.g. |Caon et al 


11991 iGraham & GuzmarJl2003l 


de Jong et alJl2004l. or in 


high-density environments (e.g. 


Gutierrez et alJl2004 IChristlein & Zabludoff 20051). In or- 



der to obtain a representative low redshift sample it is nec- 
essary to draw on a survey that is both deep and wide. 
iBlanton et al] ]2003l) fit a single S ersic function to over 
180,000 SDSS galaxies, and recently iTasca fc White! (120051) 
have performed bulge-disc decomposition on a smaller sam- 
ple of 1588 SDSS galaxies. 

In this paper we use a p ublically available b ulge-disc 
decomposition code (GIM2D; ISimard et al.l 120021) to pro- 
vide a quantitative measure of the surface brightness pro- 
files of 10095 galaxies with B M ac < 20 mag in the Mil- 
lennium Galaxy Catalogue (MGC). In Section [5] we briefly 
describe the MGC dataset analysed in this paper. Section 



|3] outlines the GIM2D fitting process. The resulting struc- 
tural catalogues are introduced in Section 13.41 and in Sec- 
tions 01 El and HJ we address the validity and repeatability 
of the structural measurements. Finally, in Section [7| we re- 
cover the fundamental empirical results describing bulges 
and discs (the Kormendy relation and the fj,o — log(/i) re- 
lation), and show that galaxy bimodality can be explained 
as a manifestation of distinct bulge and disc properties. In 
future papers we use the structural catalogues presented in 
this paper to consider the bulge and disc luminosity func- 
tions (Allen et al. in prep 2006), the bulge-disc bivariate 
brightness and size distribution (Liske et al. in prep 2006), 
and the supermassive black hole mass function (Graham et 
al. in prep 2006). Throughout this paper we assume a cos- 
mology with £l m — 0.3, £Ia = 0.7, and adopt h — _Ho/(100 
km s _1 Mpc -1 ) for ease of comparison with other results. 



2 THE MILLENNIUM GALAXY CATALOGUE 



The Millennium Galaxy Catalogue (MGC 1 ; iLiske et afl 
2003) is a deep (fj,u m =26 B mag arcsec -2 ), wide area 
(~37.5 deg 2 ) imaging and redshift survey covering a 0.5 
deg wide strip along the equatorial sky from lOh to 14h 
50'. The survey region overlaps bo th the Two-degree Field 
Galaxy Redshift Survey (2dFGRS: IColless et allboOll) . and 
the Sloan Digital Sky Survey data release 1 (SDSS-DR1; 
lAbazaiian et al J 12003) . Comparison between these surveys 
shows that the MGC (which contains 10095 galaxies with 
B < 20 mag) is deeper, more complete, more precise, 
and of higher res olution than eithe r the 2dFGRS or SDSS- 
DR1 data s e ts JCross et al] 120041: iDriver et all 120051) . In 
iDriver et al.l <l2005l) . their Figure 1 compares MGC imaging 
to DSS and SDSS imaging, and their Figure 3 demonstrates 
the extremely high redshift completeness of the MGC as a 
function of apparent magnitude, effective surface brightness, 
and colour. 

Although the input MGC imaging is in the _B-band, 
a full match to the SDSS photometric catalogues provides 
additional colour information in the ugriz bands. Full de- 
tails of the MGC, including observations, dat a reduction, 
image detection and classification, are given in lLiske et al] 
( 2003 ). The MGC redshift survey is discussed in detail in 
IDriver et al.l 12005) . With photometric precisi on of ±0.03 
mag, astrometric accuracy of ±0.08 arcsec iLiske et al] 
2003), and 96% redshift completeness to Bmgc = 20 mag 
(increasing to 99.8% for -Bmgc < 19 mag), the MGC repre- 
sents an extremely high quality and high completeness cen- 
sus of the local galaxy population. It is therefore the ideal 
data set to use for a detailed analysis of the structural com- 
position of galaxies in the local Universe, and to provide a 
low redshift anchor for higher redshift studies. 

Figure U s hows the seeing corrected half-light radii (as 
measured in lLiske et al.l l2003) for MGC galaxies as a func- 
tion of apparent B-band magnitude. Most galaxies are in- 
trinsically larger than the PSF HWHM and are therefore 
useful for structural analysis (n.b. those galaxies that are 
much smaller than the typical PSF size are eventually re- 
moved from the final analysis - see Section f5.H . 



1 http: //www. eso.org/ ~jliske/mgc 
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Figure 1. Seeing correc ted half-light rad ii (pixels) for MGC 
galaxies (as measured in iLiske et al.ll2003l) versus apparent B- 
band magnitude. The dashed line shows the typical seeing 
HWHM. The pixel scale is 0.333". 

3 DECOMPOSITION OF GALAXY PROFILES 

To perform 2D bulge-disc decomposition on t he MGC sam- 
ple w e elected to use the GIM2D package dSimard et alJ 
2002). GIM2D allows galaxies to be modelled using a single- 
component model, or a two-component bulge plus disc 
model. Where a two component model is used, the com- 
ponents are required to have a common spatial centre, but 
their luminosities are independent, allowing the calculation 
of a bulge-to-total (B/T) luminosity ratio for each galaxy. 
The intensity profile of the spheroidal (bulge) components, 
h(R), can be described using a Sersic function dSersidll968t 
iGraham fc Driverll2005l) : 

I b {R) = I e e X p(-b n [{R/R e ) 1/n - 1]), (1) 

where the effective radius, R e , encloses half the total bulge 
luminosity, and I e is the intensity at the effective radius. 
The Sersic parameter, n, determines the shape of the pro- 
file. A profile w i th n = 4 corresponds to the J? 1 / 4 model 
ide Vaucouleursl Il94gt) that was traditionally used to de- 
scribe bright elliptical galaxies. The term, b n , is defined such 
that r(2n) = 2 r y(2n,b n ), where V is the complete gamma 
function and 7 is the incomplete gamma function, ensuring 
that R e encloses half the light for all values of n . The term 
b n ca n be approximated by 1.9992n — 0.3271 llCapaccioli 
1 19891) 2 . 

For discs, the intensity profile, Id(R), can be modelled 
using an exponential function: 

2 For n < 0.36 this approximation begins to fail ICiotti fc Bertir] 
1999), leading to uncertainties of order 0.1 mag arcsec - 2 in sur- 
face brightness. U nfortunately, although m ore accurate functions 
are available (e.g. lMacArthur et al.ll200 3). this approximation is 
a fixed feature of the GIM2D code. 



I d (R)=I exp{-R/h), (2) 

where Io is the central intensity and h is the scale length 
(identical to a Sersic profile with n = 1). Although GIM2D 
permits the application of a simple opacity model for discs, 
we choose at this stage to treat discs as transparent and 
infinitely thin. We explore the effects of dust on disc opacity 
in a future paper (Allen et al., in prep). 

The ellipticity of the bulge component and the inclina- 
tion of the disc component are permitted to vary indepen- 
dently. Flattening of bulges is described by an ellipticity, 
e = 1 — b/a, where a and b are the semi-major and semi- 
minor axes of the ellipse respectively. The bulge effective 
radius, and disc scale-lengths computed by GIM2D corre- 
spond to the semi-major axes of the bulge and disc respec- 
tively. Discs can have an inclination, i, from face-on (i — 0) 
to edge-on (i = 90), measured under the assumption that 
face-on discs are circular. Bulges and discs are also permit- 
ted to have independent position angles, 0b and 8d- Finally, 
the GIM2D code allows us to recalculate the background 
levels for each galaxy before the surface brightness models 
are fitted to the galaxies. 

3.1 Galaxy Input Images and Masks 

Initia lly, GIM2D uses the SExtractor feertin fc Arnouts! 
1996) output catalogue to prepare galaxies for analysis. In 
particular, the galaxy x-y position, sky background, and 
isophotal area as determined by SExtractor are required by 
GIM2D. Postage stamp images were created for each galaxy, 
centred on the input x-y positions, with the height and 
width set to contain 10 times the SExtractor determined 
isophotal area of the galaxy. 

For each original image, SExtractor was set to produce 
a segmentation image, or mask. This consists of a pixel map 
defining background pixels to have a value of 0, and pixels 
that are considered part of objects are assigned values of 1 
to TV. A mask postage stamp image with the same dimen- 
sions as the input image was then produced for each galaxy. 
Figure |21 shows example input MGC galaxies and their as- 
sociated masks. In some cases the masks did not distinguish 
between nearby objects (in most cases two or more nearby 
objects were given the same pixel value in the mask). As we 
discuss later, this results in an erroneous output. When this 
occurred, SExtractor deblend and/or threshold parameters 
were changed to produce a correct mask (see Section f4. II for 
full details) and the input masks remade. 

3.2 Point Spread Functions 

Before accurate profile measurements can be made, it is im- 
portant to disentangle the intrinsic morphologies of galaxies 
from distortions that arise from the combined optical system 
of telescope, instrument, and atmosphere. These distortions 
vary both as a function of position and time, so their com- 
bined effects must be modelled using a point spread function 
(PSF), which is unique to each galaxy. In the fitting process, 
GIM2D convolves the model profiles with a model PSF be- 
fore comparing the model to the data. 

Initially the IRAF/DA0PH0T package was used to create 
a PSF model for each CCD frame. This was achieved using 
a minimum of 20 stars (often more), selected on the basis 
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Figure 2. GIM2D input and output for three representative example galaxies: MGC27301, MGC61361 (both shown with 
Sersic+exponential fits), and MGC55593 (single-component Sersic fit). For each galaxy the top row contains images of the galaxy 
(left), GIM2D model image (centre), and a residual image showing the difference (right). The bottom row shows the input PSF (left, 
note the expanded scale), the SExtractor segmentation image or mask (centre), and a plot of the best-fitting profile found by G1M2D 
(right). 



© 0000 RAS, MNRAS 000, 000-000 



MGC: Bulge/Disc Decomposition of 10095 Galaxies 5 



of their SExtractor stellarity parameter being > 0.8, and 
their magnitude lying in the range 17 < Bmgc < 19.5 mag. 
The stars were chosen to provide good spatial coverage of 
the CCD, to not be blended, to have no bright neighbours 
or saturated pixels, and to be away from CCD defects. The 
DAQPHOT routine PSF was then used to fit a model PSF for 
each CCD, based on the surface brightness profiles of the 
stars. 

The PSF for each image was modelled using an analytic 
function, and a lookup table of residuals to account for the 
variation of the PSF across the CCD frame. The PENNY2 
function, which consists of an elliptical Gaussian core with 
Lorentzian wings (both of which have independent position 
angles) was used to model the PSF, with the fitting radius 
of the function (i.e. the radius which encloses the region of 
pixels used in fitting the analytic function) defined to be 
equal to the median seeing for the frame. Where the me- 
dian seeing is the median FWHM of simple Gaussian fits to 
the profiles of stars. The diameter of the region used for fit- 
ting the PENNY2 function is therefore twice the estimated 
FWHM. The PENNY2 function consists of five free param- 
eters, P n (where n = 1...5), and a normalisation factor, A, 
and can be expressed as: 

PENNY2 = A x {^y^j + Ps exp(-0.693e)) , (3) 

where 

x 2 y 2 

2 = pS + h +xyPs, (4) 
and 

x 2 y 2 

e =-pI + -52+ x V p i- ( 5 ) 

The resulting PSF image is the best-fitting function to the 
light profile of the stars in the frame weighted by their signal- 
to-noise. 

To check the quality of the PSF, the ALLSTAR routine 
was used to reconstruct and subtract the PSF from the stel- 
lar images. When the PSF is good, the stars are cleanly re- 
moved from images leaving only sky noise. In most cases the 
PSF provided a good model for the stars (x?ed ~ 1) where 
Xred is the reduced \ 2 f° r the fit), although for a small num- 
ber of frames the Xred values were poor and additional stars 
were rejected, and the PSF regenerated until an acceptable 
Xrod was obtained. 

Using the model PSFs made for each CCD image, with 
the associated look-up tables, the DA0PH0T routine SEEPSF 
was employed to create an image of an artificial star corre- 
sponding to the exact location of every MGC galaxy. GIM2D 
then convolves this image with the model fits before com- 
parison with the input data. Example PSF stars are shown 
in Figure [5] (note the expanded scale for the PSFs). 

3.3 The GIM2D Models 

Three different models were fit to each MGC galaxy. 

• One component: fitted with a Sersic function. 

• Two components: de Vaucouleurs bulge + exponential 
disc. 

• Two components: Sersic bulge + exponential disc. 



A full description of t he GIM2D fitting algorithm is 
given in lSimard et al.H2002l) . but a brief outline follows. The 
first stage in the fit is to further refine the background sub- 
traction. GIM2D uses those pixels flagged in the SExtractor 
segmentation image as 'background' to estimate the mean 
value of the background. A 5 pixel buffer is placed around 
each object in the segmentation image to exclude possible 
object pixels from isophotes lower than the threshold that 
could potentially bias the background estimation. The new 
background value is then fixed, and 2D surface brightness 
fits are applied to the galaxies. 

Table Q shows the 8-12 parameters that are fit by 
GIM2D and the the upper and lower limits permitted for 
each of the fits, most of which are based on GIM2D de- 
faults. Since some galaxies are known to have n > 4.0 (e.g. 
ICaon et al 1ll993t iGraham et"aflll996h. we permit th e sersic 
index to be as high as 12 fc.f. lTasca fc Whitdl2005l who fix 
Mupper = 4). Since we have assumed an infinitely thin disc, 
inclinations are given an upper limit of 85 degrees. Bulge 
ellipticities are permitted to have a maximum value of 0.7. 
However, in the Sersic-only catalogue the fit may best cor- 
respond to either a bulge or a disc. Therefore, in this cat- 
alogue, those objects with a GIM2D upper error limit of 
0.7 for the ellipticity were refit using a revised upper limit 
of 0.92 (corresponding to an inclination of 85 degrees), and 
both versions of the fit stored. 

GIM2D uses the pixels flagged by the segmentation im- 
age as 'object' to measure image moments in order to pro- 
duce initial estimates of the total flux, size, and position 
angle for each galaxy. Maximum values for these parame- 
ters are then set at twice the best estimates from the im- 
age moments. Using the limits applied in Table and the 
maxima derived from the image moments, N models are 
created coarsely sampling the permitted parameter space. 
For this analysis we set N = 400. An 'initial condition 
filter' is then used to search through these models to find 
the best fit. The values associated with this model are then 
used as the starting point for the final fit. GIM2D then uses 
the Metropolis algorithm to search for the maximum like- 
lihood model (essentially a pixel-to-pixel \ 2 minimisation) 
that best describes each galaxy. GIM2D also explores the 
parameter space around the best fit to compute 68% confi- 
dence intervals for each parameter. 

3.4 MGC Structural Catalogues 

Publically available MGC catalogues have already been re- 
leased containing non-GIM2D parameters, such as photom- 
etry, imaging, estimates of galaxy half-light radii, mean 
effective surface b rightnesses, and a match to SDSS-DR1 
jLiske et al]l2003l) . These were supplemented w ith redshifts 
and rest-frame colours in iDriver et alJ i2005T) . The best- 
fitting values for the twelve parameters computed by GIM2D 
(see Table [TJ have been combined with the MGC-BRIGHT 
catalogue (i.e. 10095 galaxies with B < 20 mag) to produce 
three structural catalogues: one for each of the three models 
fitted here. For each fit, GIM2D also measures the reduced 
X 2 (as computed between the PSF-convolved best-fitting 
model and the input data) as a measure of goodness-of- 
fit. Furthermore, using the topology of parameter space ex- 
plored during the fitting process, GIM2D comp utes 68% con- 
fidence limits for each of the parameters (see lSimard et ail 
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Figure 3. Residual offsets in the location of the object centre as determined by SExtractor and GIM2D (in x-y pixels), are shown for the 
de Vaucouleurs + exponential catalogue (left), the Sersic + exponential catalogue (centre), and the single-component Sersic catalogue 
(right). The dashed lines correspond to a difference of 3 pixels (= 1 arcsec) in either x or y. 



Parameter 




Lower 


Upper 


Total Flux 


F (counts) 


0.0 




Bulge/Total Flux Ratio 


B/T 


0.0 


1.0 


Bulge Effective Radius 


R e (pixels) 


0.0 




Sersic Index 


n 


0.2 


12.0 


Bulge Ellipticity 


e 


0.0 


0.7 


Bulge Position Angle 


B\y (degrees) 


-360.0 


+360.0 


Disc Scale Length 


h (pixels) 


0.0 




Disc Inclination 


i (degrees) 


0.0 


85.0 


Disc Position Angle 


6d (degrees) 


-360.0 


+360.0 


x —Position 


A x (pixels) 






y —Position 


A y (pixels) 






Background level 


A b (counts) 







3.4-1 Other Parameters 

After profile fitting has taken place and output and resid- 
ual images made, GIM2D also computes several measures 
of con centration and asym metry which are described in de- 
tail in lSimar d et alj (|2002|). These include the C — A system 
(based on TATr^Jianiet^ilJl l99(J) and, using residual images, 
Rt and Ra JSchade et al.lll995T). Two other measures, A z , 
and D z , are also defined bv ISimard et al] i2002t) as part 
of the GIM2D package. Although these parameters are not 
discussed further in this paper, the C — A parameters are 
made available in the final MGC structural catalogues, and 
the other asymmetry measures are available on request. 



Table 1. The 12 parameters that can be fitted by GIM2D, and 
the upper and lower limits permitted in the fits. No entry indi- 
cates that no hard limits are applied. For the Sersic-only (one 
component) fit, the B/T lower limit is constrained to 1.0, so only 
the Sersic component is fit. When a de Vaucouleurs (R 1 / 4 ) fit is 
applied, the upper and lower limits of the Sersic index are held 
fixed at ra = 4.0. 



2002). Although these errors are included in the catalogues, 
they only account for uncertainties in the formal fits. Un- 
certainties and errors due to the sky background and PSF, 
which can often dominate, are not considered in these esti- 
mates. This is discussed further in Section E] 

For each galaxy, in addition to the catalogue of struc- 
tural parameters, output PSF convolved model images, and 
residual images are also produced by GIM2D. Figure 
shows example output and residual images for three galax- 
ies: MGC27301, MGC61361, and MGC55593, along with the 
raw images, masks and PSFs, and a plot of the best fitting 
profiles. Further details of our final catalogue parameters are 
presented in Appendix A 3 . 



Structural catalogues and images are made available at the 
MGC website http: //www. eso.org/ ~jliske/mgc 



4 QUALITY CONTROL 

4.1 Comparison with Independent Measurements 

As an initial test on the accuracy of the GIM2D output, sim- 
ple global observables such as magnitudes, half-light radii, 
and x — y centroid positions can be co mpared to th e ir val - 
ues measured using other means in the iLiske et"ai1 (120031) 
MGC catalogues. In most cases the agreement is excellent, 
although there are a number of outliers, which can be a 
useful part of the quality control process as we discuss in 
Section IO 

Figure [3] shows the residuals between the SExtractor 
and GIM2D centroids. The majority of objects (88.0%) lie 
within 1 pixel of the SExtractor central position, and 99.1% 
have GIM2D central positions within 1 arcsec (3 pixels) 
of the SExtractor position for the Sersic+exponential cat- 
alogues, and similar accuracy is found for the other two 
models used in this paper. In most cases an outlying ob- 
ject is simply a sign of an irregular galaxy with an offset or 
ill-defined core, although it can, as we discuss in the next 
section, indicate problems with the masks (SExtractor seg- 
mentation images), or the GIM2D output. 

Figure 2] shows the difference between the ILiske et alJ 
( 2003) and GIM2D magnitudes as a function of MGC mag- 
nitude, and the best-fitting relations between these parame- 
ters along with residual histograms, and the mean and stan- 
dard deviation (with iterative 3a clipping). For all three 
structural catalogues the GIM2D magnitudes are brighter 
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Figure 4. Differences between MGC galaxy magnitudes (de- 
termined by SExtractor) and GIM2D magnitudes, as a func- 
tion of MGC magnitude. The solid line shows the least squares 
fit to the data, and the dashed lines define the 1 — a enve- 
lope. There is no apparent trend with magnitude. The GIM2D 
magnitudes are slightly brighter in all three catalogues: de Vau- 
couleurs+exponential (top), Sersic+exponcntial (middle), and 
Sersic (bottom). The panels on the right show histograms of the 
residuals. The mean and standard deviation (with 3ct clipping) of 
these distributions arc shown in each case. 



by 0.07 — 0.09 mag, with a 0.05 — 0.07 mag standard devi- 
ation. There appears to be no significant trend with mag- 
nitude. Systematically brighter GIM2D magnitudes are not 
unexpected since profiles are integrated to R = oo, and it is 
known that flux can be missed when usin g aperture or Kron 
magnitudes (see lGraham fc D river 2005). There are a small 
percentage of objects with large residuals that we discuss in 
the next Section. 

Finally, GIM2D also computes an approximate measure 
of the total galaxy half-light radius, which assumes that the 
bulge and disc have the same position angle (although the 
position angles are independent in the fits). In Figure the 
GIM2D galaxy half-light radii are com pared to th e seein g 
corrected half-light radii measured in iLiske et alJ {2003). 
Again, there is good agreement but with some scatter. It 
is also notable that at small radi i (< 1 .2") the GIM2D radii 
are smaller than the lLiske et alJ J2003) radii. This is unsur- 
prising as it is in the regime where the PSF is of comparable 
angular size to the galaxies. GIM2D measurements will be 
dominated by errors due to the PSF correction and the orig- 
inal MGC measurements will be su bject to the assump tions 
of the seeing correction applied bv ILiske et alJ i2003l) . For 
larger, well resolved galaxies the GIM2D half-light radii are 
slightly larger, as expected, given their brighter magnitudes. 



4.2 Reanalysis of Outliers 

Analysis of the input and GIM2D output images of outlying 
objects in Figures [3] and 2] reveals that the difference in pa- 



rameters between ILiske et alJ (l2003h and GIM2D measure- 
ments can have a number of sources. In some cases galax- 
ies are irregular, with more than one possible 'centre', or 
have a very asymmetric surface brightness distribution. The 
GIM2D models can be a poor fit to such objects, and po- 
sition and magnitudes consequently have significant differ- 
ences. In addition, many objects have nearby neighbours, in 
which case SExtractor has applied a deblending algorithm to 
attempt to separate the flux from the two objects. SExtrac- 
tor creates a mask image with overlap pixels simply desig- 
nated to one object or the other. However, when SExtractor 
computes the magnitudes for these objects, the flux in the 
pixels lying in overlap regions is shared between the two ob- 
jects. No such deblending is performed by GIM2D, which 
simply uses the mask, and all the flux is assumed to belong 
to one object or the other. Therefore some GIM2D inputs 
can be contaminated by light from nearby neighbours, and 
the derived magnitudes will differ from those determined by 
SExtractor. 

Another cause of large parameter differences are erro- 
neous masks or segmentation images. SExtractor, like most 
automated detection algorithms, is typically set to run with 
'optimal' tuning (i.e. the detection parameters are set so 
that reasonable output is produced for as many objects as 
possible) . Nevertheless, a fraction of objects at the extremes 
will still be incorrectly recovered. Perhaps the best examples 
are the SExtractor deblending parameters DEBELND.NTHRESH 
and DEBELND_MINC0NT which govern the number of deblend- 
ing thresholds, and the minimum contrast used when SEx- 
tractor decides whether two nearby local maxima in sur- 
face brightness are significant enough to be considered sep- 
arate objects, or just structure within a single object. At 
one extreme, if the deblending is too strong, objects such 
as spiral discs or irregular galaxies can be broken up into 
several smaller objects. At the other extreme, two nearby, 
but clearly distinct objects are considered one, and assigned 
the same arbitrary pixel value in the mask (this is especially 
common when a small object lies close to a much larger 
one). To confuse matters further, in a small number of cases, 
SExtractor has failed to assign different pixel values to two 
nearby objects in the mask, even if two objects appear cor- 
rectly deblended in the catalogue. These errors can be fixed 
by changing the deblend parameters to be more or less sen- 
sitive depending on the problem, and rerunning SExtrac- 
tor. However it is not possible, even using the most optimal 
setup, to produce output that is satisfactory in every case. 
This is a problem when trying to automate measurements 
for increasingly large data sets as it becomes more and more 
difficult to check all objects by eye (and defeats the purpose 
of trying to automate the process!). Either some fraction of 
wrong detections has to be tolerated, or a way needs to be 
found where the failures can be flagged and corrected. 

For the MGC, a ll 10095 objects were checked by eye 
in ILiske et alJ <l2003l) . a nd therefore , as al ready discussed, 
differences between the ILiske et alJ B200a) catalogues and 
GIM2D output can be indicative of bad masks. If a galaxy 
has been over-deblended, then the GIM2D magnitude will 
be an under-estimate. When two galaxies are treated as one, 
the centroids will be grossly different, and the GIM2D mag- 
nitudes over-estimated. All objects with Ax or Ay greater 
than 3 pixels (1 arcsec) or a magnitude difference | AB| > 0.2 
mag, were selected and inputs, outputs, and masks examined 
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Figure 5. Comparison between galaxy half-light radii measured by GIM2D, and seeing corrected half-light radii computed in lLiske et alj 
(2003), for the de Vaucouleurs+exponential catalogue (left), the Sersic+exponential catalogue (centre), and the single component Sersic 
catalogue (right). There is good agreement except at small radii where the PSF dominates. 



by eye. It was found that 222 masks needed to be corrected 
by hand. 

An example is shown in the top row of Figure |S| with 
the original input image, bad mask (which assigns the same 
pixel value to two distinct objects), and the clearly erroneous 
GIM2D fit on the right. A corrected mask, which is produced 
by making small adjustments to SExtractor deblending and 
threshold parameters, is shown in the bottom panel, along 
with the (now corrected) GIM2D output. 



5 INTERPRETING GIM2D OUTPUT 

One of the problems with automatic bulge-disc decomposi- 
tion routines is that in addition to perfectly good fits, they 
automatically generate a lot of rubbish. In this Section we 
describe our best efforts to clean up this situation. 

As a starting point we first trim the catalogue and ex- 
amine only those systems within the redshift range 0.013 < 
z < 0.18, and impos e the galaxy s i ze and surface brightness 
limits introduced in iDriver et alj <l2005l). This results in a 
trimmed catalogue of 7750 galaxies (see lDriver et al.ll2005(l 

The interpretation of single-component Sersic fits is 
fairly trivial, with the resulting light profile assumed to be 
the best-fit to the light from the entire galaxy. However, 
when two components are modelled, the two functions are 
typically taken to correspond to a bulge and a disc. It is 
necessary to verify whether or not this is the case, ensuring 
that the two-component model really does correspond to 
two distinct structural components within a galaxy. This is 
especially important in the case where a Sersic+exponential 
model is used because this fit has the most degrees of free- 
dom. 

No restrictions on the relative sizes of the 'bulge' and 
'disc' components were applied when using GIM2D. The two 
components can sometimes be inverted or used to fit other 
structural features, especially in irregular galaxies. The com- 
ponent light profiles may also cross once, twice or not at all. 
To ensure the correct interpretation of the GIM2D fits (and 
whether components should be interpreted as bulges, discs, 
or something else) we classify all profiles into one of eight dif- 
ferent types: six are composed of two components, and two 
have only single-component profiles. The two-component 



profile types are plotted in Figure|7|for six example galaxies 
from the Sersic + exponential catalogue. The profile types 
can be summarised as follows: 

• Type 1 'Classic' profile. The Sersic profile dominates 
at the centre, while the exponential profile dominates the 
flux at large radii. At a surface brightness brighter that 26 
B mag arcsec -2 the profiles cross only once (i.e. they have 
the same surface brightness at only one radius). 

• Type 2 Disc dominated system. The exponential pro- 
file dominates at all radii, with a small, central Sersic profile. 
The profiles never cross, and B/T is always < 0.5. 

• Type 3 Sersic profile dominates at small and large 
radii, but an exponential profile dominates at intermediate 
radii. The profiles cross twice. B/T is always > 0.5, and the 
Sersic index, n is typically > 1.5. Objects fitted by these 
profiles are typically single-component Elliptical galaxies. 

• Type 4 Inverted profile. The exponential profile domi- 
nates at the centre, and the Sersic profile dominates at large 
radii. The profiles cross only once, and n is always small. The 
correct interpretation of the profile is problematic, and could 
be a signature of possible disc truncation, poor background 
subtraction, or the profile could correspond to irregular or 
dwarf systems. 

• Type 5 Bulge or disc dominated. The Sersic profile 
dominates at all radii, with a weak underlying exponential 
component. The profiles never cross, and B/T > 0.5. There 
is a bimodal distribution in n for this type (see Section f5. Ill . 

• Type 6 Disc with perturbation. The exponential pro- 
file dominates at small and large radii, and the Sersic profile 
briefly dominates at intermediate to large radii. The Profiles 
cross twice. The Sersic component models either a disc per- 
turbation (such as a spiral arm), or features in an Irregular 
galaxy. 

• Type 7 Pure disc: an exponential-only profile with 
B/T = 0. 

• Type 8 Pure bulge, disc, or other: Sersic profile only 
with B/T = 1. A bimodal distribution in n reveals discs 
with n ~ 1, and bulges with n ~ 4 for our sample. 

All 8 profile types occur in the Sersic+exponential 
catalogue, types 1,2,3,5,7 and 8 occur in the de Vau- 
couleurs+exponential catalogue, and the Sersic-only cata- 
logue only contains type 8 profiles by definition. In both 
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Figure 6. An example of bad and corrected mask images. The top left panel shows the input image for MGC64481 (the central object), 
and the original mask (top centre). The top right image is the original GIM2 D output image, wh ich can be initially identified due to its 
large centroid offset, and magnitude difference between the GIM2D output and lLiske et al.H2003ly The bottom panels show the corrected 
mask and output images. 
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Table 2. The fraction of galaxies corresponding to each profile 
type (see Figure 171. 



the Sersic+exponential and the de Vaucouleurs+exponential 
catalogues the majority of galaxies are of type 1. In the de 
Vaucouleurs+exponential catalogue ~ 1/3 of galaxies are 
best modelled by an exponential-only fit (type 7), three 
times as many as in the Sersic+exponential catalogue. Table 
[5] summarises the fraction of each type that is found in the 
three different catalogues. 

5.1 The Logical Filter 

The use of raw GIM2D profiles from the two-component 
models would clearly be inappropriate given the distinct and 
often physically meaningless profiles that are fitted in some 
cases (e.g., types 4 or 6). Moreover, many systems do not 
actually exhibit two components, and in these cases it is not 



surprising that a two-component fit produces erroneous out- 
put. A non-negligible fraction of incorrect output is typical 
for most automated algorithms. A quick in spection of the 
bulge -disc decomposition s from BUDDA jde Souza et alJ 
l2004|) . and lAguerri et all ll2005t) reveal that problems with 
automatic decompositions are pandemic, with many fits hav- 
ing unrealistic bulge-to-disc size ratios, disc components 
clearly falling faint of the outer discs in spiral galaxies, and 
bulge flux dominating at large radii in late-type spiral galax- 
ies, etc. While automatic codes find the optimal mathemat- 
ical solution, perturbations in the real data from the fitted 
models often result in the mathematical solution being an 
unphysical one. 

However, it is possible to ask whether one can iden- 
tify a suitable set of 'rules' to apply to the GIM2D cat- 
alogues to produce a meaningful, final 'filtered' catalogue. 
This catalogue would consist of bulge-disc decompositions 
where appropriate (using the Sersic+exponential fits), and 
single-component Sersic fits otherwise. The interpretation 
of the single-component Sersic fits (as a bulge or a disc) 
would depend on colour and/or Sersic index. This same ap- 
proach could be followed using the de Vaucouleurs catalogue 
for bulge-disc decomposition but we choose to consider only 
the more meaningful Sersic+exponential catalogues here. 

For each fit GIM2D computes the reduced-x 2 as 
a measure of the goodness-of-fit for the output model. 
In most cases both the Sersic+exponential fits and 
the single-component Sersic catalogues have reasonable 
reduced-x 2 values (~ 1), although the Sersic+exponential 
reduced-x 2 values are generally smaller. Although the 
Sersic+exponential fits may have a reasonable reduced-x 2 , 



© 0000 RAS, MNRAS 000, 000-000 



10 Allen et al. 




Figure 7. Examples surface brightness profiles for the 6 distinct types of GIM2D fit that we use to classify galaxies from the Sersic + 
exponential catalogue (see text). In each panel the bold solid line represents the total profile, the Sersic component is represented by a 
dotted curve, and the dashed straight line represents the exponential component. The locations of the half-light radii for each component: 
R e , and 1.678ft, are also shown. Table |3 summarises the frequency with which each type occurs. 



this doesn't indicate whether or not the fit is appropriate 
(i.e. whether the Sersic function is really corresponding to a 
bulge and the exponential function to a disc) . Therefore the 
approach followed here involves analysis of the global prop- 
erties of each type (including colours, Sersic indices, and 
the sizes of components), along with visual examination of 
examples of each of the eight types. 

Objects have a range of B /T values (Figure |H] shows 
the B/T distribution for our type 1 galaxies), although we 
note that the frequency of high B/T values is small (in- 
deed, only 62/7750 objects are type 8 with B/T = 1). This 
perhaps reflects the fact that the addition of a weak ex- 
ponential model to a light profile dominated by a Sersic 
model can make little differe nce to the overall total pro- 
file (see e.g. lSaglia et al.l [l997). The exponential component 
may represent a weak disc, an isophotal twist, a PSF error, 
or some perturbation in the profile. This degeneracy means 
that kinematic measurements are required to test for the 
presence of a low- luminosity, rotationally supported disc. 

Figure [5] shows the distribution of disc inclinations for 
spheroid dominated systems (type 1 and type 5 galaxies with 
n > 1.5) for B/T < 0.8 and B/T >= 0.8. For objects with 
B/T < 0.8, the distribution is relatively flat, as expected for 
a random sample of discs with random orientations in space. 
However, for B/T >= 0.8, 'discs' with inclinations towards 
face-on dominate, especially for those objects that have been 
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Figure 8. The distribution in B/T for (classic) type 1 galaxies. 
The lack of objects with very high B/T is notable. 
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Figure 9. Distribution in cos(i) for the 'disc' component of galax- 
ies with types 1 and 5 with a bulge component having n > 1.5. 
The top panel shows a roughly flat distribution for those galax- 
ies with B/T < 0.8, whereas the lower panel, for galaxies with 
B/T > 0.8, indicates a suspicious preference for face-on exponen- 
tial profiles. 



visually classified as E/S0. It is therefore highly improbable 
these faint discs are all real. Therefore, we choose to replace 
the Sersic + exponential fits for type 1 fc 5 objects with 
B/T > 0.8 with single Sers ic profiles (c.f. lTruiillo fc Aguerril 
l2004l:lGutierrez et all2004 who replace objects with B/T > 
0.6 with single Sersic p rofiles). 

iDriver et alJ i2006T) have demonstrated the existence of 
a bimodality in the colour-rog(n) plane for luminous galax- 
ies (when modelled with a single Sersic function); one can 
clearly identify a high-n red peak, and a low-n blue peak. 
The colours referred to here are the (u — r) COTe PSF colours 
from SDSS photometry (see lDriver et al.ll2006t for more de- 
tail of our use of PSF colours). Furthermore, they suggest 
that the two peaks can be identified with bulges and discs, 
respectively. This differs in a subtle but important way from 
the usual interpretation of the colour bimodality, which is 
credited to early- and late-type galaxies rather than bulge 
and disc stellar systems. Figure [TU1 shows the colour-log(n) 
plane for the Sersic (which one might initially interpret as a 
'bulge') component for each of the six two-component types 
we identify (see Figure |7J. The top-left panel of Figure [Til 
shows the same cumulat ive raw distribution for all types to- 
gether. This differs from lDriver et aD i2006T) in that we have 
now separated the bulge and disc components, as Driver et 
al. argued should be done. It is clear from Figure [TU1 that 
most galaxies are 'classical' type 1 profiles, and of these, 
the majority of the 'bulges' are associated with a high-n red 
peak. However, for most galaxies of types 2-6, and for a 
small but significant fraction of type 1 galaxies, the bulges 
do not correspond to the high-n red peak. To investigate this 
situation further we now consider each of the profile types 
in turn. 

The majority of the two-component fits correspond to 



Figure 12. The distribution of bulge Sersic index, n, for type 1 
galaxies that are bulge dominated (i.e. B/T > 0.5). 



a 'classic' type 1 profile. Analysis of the bulge properties 
of these galaxies suggests the presence of at least two pop- 
ulations in the majority of cases. The top left panel of Fig- 
ure 1101 shows that most objects lie in a clear peak around 
n ~ 3, and PSF colour (it — r) colo ~ 2.5. These bulges 
appear to belong to plausible bulge-disc systems, and are 
accepted by us as correct decompositions, and therefore 
progress unchanged to the final filtered catalogue. Addition- 
ally, there is a tail of objects and a possible second peak 
around n ~ 0.5, with a wider range in colours, albeit gen- 
erally blue [(it — r) corc < 2.0]. This could be indicative of 
a population of pseudo-bulges or bars. Such a population 
might be expected in disc dominated systems {B/T < 0.5), 
and these objects are accepted as having good decompo- 
sitions and genuine bulges. An additional separation into 
bulges and pseudo-bulges/bars can be made at a later stage 
if required. However, pseudo-bulges are unlikely to be found 
in the high B/T regime. Figure WH shows the distribution 
in Sersic index for Sersic dominated (B/T > 0.5) type 1 
galaxies. Based on this distribution we accept objects with 
n > 1.5 as having a good bulge-disc decomposition and they 
are included as such in the filtered catalogue, whilst type 1 
objects with n < 1.5 are replaced with single-component 
Sersic fits. These objects appear to be either pure disc sys- 
tems or dwarf ellipticals. 

Type 2 galaxies are clearly disc-dominated objects 
(they all have B/T < 0.5), with a possible small bulge. The 
addition of a weak Sersic function to an exponential func- 
tion can be enough to perturb the total profile significantly 
(see the type 2 example in Figure[7J. However, the modelled 
perturbation to the total profile is not always at the cen- 
tre, and is sometimes at much larger radii (analogous to the 
type 6 profile shown in Figure |7| but without the profiles 
crossing). In these cases it would be incorrect to interpret 
the Sersic function as a bulge, and we therefore impose the 
additional restriction for type 2 objects that R e < 0.5 x h, 
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Figure 10. The distribution of galaxy bulges (pre-filtering) in the (u — r) cor c — log(n) plane for each of the six profile types we identify 
in Section l5l 

ensuring that only Sersic functions modelling a central fea- Many type 4 profiles can be considered true inversions, 

ture are considered as bulge-disc systems. We allow these with the exponential function fitting a central feature, and 

objects to pass through to the filtered catalogue. Based on the Sersic function fitting the outer disc. This can occur for a 

the colour-log(n) distribution in Figure tTTjl the bulges of the number of reasons. A PSF error can mean that the exponen- 

majority of these galaxies do not correspond to the classical tial is used to fit a sharp 'spike' at the centre of the profile, 

red, high-n peak, and may be better interpreted as pseudo- Such spikes are always much smaller than the HWHM of 

bulges. A small fraction of these objects are red, but they are the PSF. In other cases, the galaxy appears to be an irreg- 

all highly inclined systems, and thus likely red due to dust. ular disc with active star formation. The size and shape of 

Those objects with R e > 0.5 x h (i.e. the Sersic function the two components can be comparable, especially for the 

does not correspond to a central perturbation) are replaced redder galaxies, when the modelled galaxy is often a dwarf, 

with single-component Sersic fits, and are interpreted as In all cases the bulge-disc interpretation would be erroneous 

disc-only systems. and a single Sersic fit is preferred in the filtered catalogue. 

Only a small fraction of galaxies are classified as type However a fraction of these galaxies appear to be genuine 
3. These fits are erroneous and are sometimes caused by bulge-disc systems with severe disc truncation. This seems 
nearby neighbours (which GIM2D would be unaware of) to be a real phenomenon, and is likely to be an increasing 
contributing extra flux at large radii. They are generally problem with future, deeper data sets. Similarly, at higher 
red, bulge-dominated systems, with high-n and we choose redshifts, discs are s een to be described wit h a Sersic func- 
to replace them with a Sersic-only fit, which will typically tion requiring n < 1 llTamm fc TenieJ2006l. Unfortunately, 
be later interpreted as a bulge-only system. This class does modelling disc truncation is beyond the scope of GIM2D, 
contain a fraction (< 1%) of genuine bulge-disc systems that and it is not possible to fit the 2D images accurately. There- 
will be missed. fore, we choose to replace these fits with single-component 
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Figure 11. Left column: the distribution of galaxy 'bulges' in the (u — r) colo — log(n) plane for: (top left) raw galaxy bulges (pre- 
filtering) of all types (see Figurc llOl . (middle left) bulges from two-component Sersic + exponential model that pass through the logical 
filter unchanged, (bottom left) replacements from the Sersic-only catalogue that are interpreted as bulges. The right panels show the 
distribution of bulges from the final filtered catalogue (reliable Sersic + exponential bulges, plus replacement Sersic-only bulges) in 3D 
(top), and in 2D (bottom). 



Sersic functions, and accept that a small fraction (<5%) 
of bulge-disc systems may be missed, or rather, not decom- 
posed. 

The type 5 galaxies all have B/T > 0.6, and are there- 
fore dominated by the Sersic component. The addition of an 
exponential function (that does not cross) makes little dif- 
ference to the total profile. The distribution in Sersic index 
and colour is clearly bimodal (figure HITl , and there is a corre- 
lation between (two-component) Sersic-'bulge' parameters, 
and the Sersic-only parameters. These galaxies tend to be 
blue, low-n disc-only systems, or red, high-n bulge-only sys- 
tems (i.e. spheroids). They are best fitted by a single Sersic 
component, and these fits are therefore used in the filtered 
catalogue. 

The fits to type 6 galaxies are clearly incorrect if in- 
terpreted as bulge-disc decompositions. From Figure [THl the 
'bulge' components have very low n, and blue (it— r) colours. 
These galaxies also have low Sersic index (n < 1) from the 



Sersic-only fits. They appear to be either dwarf galaxies or 
very flat (n < 1) discs with active star formation and/or 
other irregularities. In all cases a single Sersic fit is adopted 
as more appropriate. 

Finally, those objects that only have one component 
from the Sersic + exponential model fits (type 7 and type 
8) are replaced with single-component Sersic-only fits. 

In summary, all profiles were replaced with Sersic-only 
fits except for types 1 and 2, which generally keep the two- 
component bulge-disc decomposition. The exceptions are 
when (i) B/T > 0.8, or (ii) a low-n (< 1.5) Sersic com- 
ponent dominates (B/T > 0.5) for the type 1 galaxies, and 
(iii) when R e > 0.5 x h (i.e. the Sersic function models a 
non-central perturbation) for the type 2 galaxies. The logi- 
cal filter is summarised in Figure [Til After the filtering pro- 
cess, the full catalogue contains 3454 two-component objects 
where the Sersic component is treated as a bulge, and the 
exponential component is treated as a disc. The remaining 
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Figure 13. The distribution of those galaxies replaced with 
Sersic-only fits by the logical filter in the (u — r) coro — log(n) 
plane. We choose n = 1.5 as the division between the two popu- 
lations. 

galaxies in the catalogue (4296 objects) are fitted with a sin- 
gle Sersic component. The distribution of these Sersic-only 
objects in the colour-n plane is shown in Figure Hoi where 
the bimodal nature is evident. Based on this distribution we 
choose n — 1.5 as the division between the two populations. 
We classify 825 galaxies with n >— 1.5 as bulge-only ob- 
jects, and 3471 galaxies with n < 1.5 as discs. As already 
discussed, these replacements will include a small fraction 
(< 6%) of type 3 or type 4 galaxies that are genuine bulge- 
disc systems, that are misinterpreted here as bulge-only or 
disc-only. 

As a useful verification of this bulge- or disc-only inter- 
pretation we also looked at the above Sersic-only replace- 
ments in the ii 1 ' 4 + exponential catalogue. For the replace- 
ment 'disc-only' objects 80.1% of galaxies had a best-fitting 
R 1/4 + exponential model with B/T = (i.e. GIM2D found 
a single exponential disc the most optimal model). The re- 
mainder of the replacement discs were all in disc-dominated 
systems with B/T < 0.5 (half of them with B/T < 0.1). 
For the replacement 'bulge-only' objects, the R 1 ^ + expo- 
nential model produced mostly high B/T output, with a 
median of B/T ~ 0.8, suggesting that this population is 
indeed dominated by E/SO galaxies. 



6 REPEAT OBSERVATIONS AND 
PARAMETER ACCURACY 

The geometry of the MGC is such that each pointing of the 
INT Wide Field Camera overlaps by approxim ately 0.027 
deg 2 with the previous one (see Figure 2 in iLiske et alJ 
2003), building up a mosaiced strip across the sky. Ob- 
jects that lie in these overlap regions are therefore observed 
twice. Although the MGC catalogues only utilise one of these 
observations, the imaging data exists and SExtractor cata- 



logues, and GIM2D input images were prepared for all 702 
duplicate galaxies. Here we remove those objects with incor- 
rect masks (see Section [4.211 . reducing the number of twice- 
observed comparison objects to 682. 

In each case, repeat observations lie on a different CCD 
of the WFC, with different PSFs, airmasses, sky bright- 
nesses, noise, and seeing conditions. Often, the observations 
were carried out on a different night or even as part of dif- 
ferent observing runs several months or years apart. They 
therefore provide an excellent test of the repeatability of 
GIM2D fits for a diverse and representative sample of galax- 
ies under different conditions. 

In all three catalogues acceptable repeatability (most 
measurements agree better that 20%) is obtained if a 
bulge/disc component size limit is imposed such that the 
half-light radius is larger than the half-width half maximum 
of the seeing disc (half-light radius > 0.5 x seeing). Bulges are 
compared when i? e > 0.5F, and discs are compared when 
1.678/i > 0.5r, where F is the average seeing for the image 
the galaxy lies in. Below this threshold the differences be- 
tween repeat measurements begin to increase. In addition 
we also find that the repeatability is poor for less luminous 
components, and we therefore impose a M = — 17 B mag 
limit on the component magnitudes. Note that the uncer- 
tainties discussed here are not included in the final structural 
catalogues which only contain the 68% confidence intervals 
computed by GIM2D (see Section l3~4l . 

Profiles from the repeat observations have been inter- 
preted according to the rules described in Section [5. II and 
replacements made where necessary. In most cases the profile 
type is the same for each pair of observations of an object, 
but where the profile types differ, components interpreted as 
bulges or discs are directly compared (e.g. even if a galaxy 
has a type 2 profile in one measurement and a type 1 profile 
in the other, the bulge and disc components are still com- 
pared). When one measurement of an object has resulted 
in a single-component Sersic profile and the other measure- 
ment uses a two-component profile (e.g. type 1), the Sersic- 
only profile is interpreted as either a bulge or a disc (see 
Section f5. It . and compared with the corresponding compo- 
nent in the other fit. Type discrepancies such as this occur 
about 12% of the time. For each parameter we compute the 
mean residual and the 3— a clipped standard deviation. In 
all cases the size of residuals correlates with the apparent 
half-light radii of the components. 

Figure ITS! shows a comparison between bulge and disc 
parameters from repeat observations for the filtered cat- 
alogue using two-component Sersic +exponential fits (or 
Sersic replacements). In all cases the mean is close to zero. 
We find bulge R e differences have a standard deviation of 
12.2%, the Sersic index, n, has a standard deviation of 
Alog(n) = 13.2%. Ellipticity differences have a standard 
deviation of 4.1%, and the bulge magnitudes, Mb u i go , have 
a standard deviation of 0.10 B mag. For discs a larger sam- 
ple can be used, as more discs are larger than the 0.5r cut. 
Scale lengths, h, have a standard deviation of 6.6%, and disc 
inclinations, cos(z), have a standard deviation of 4.7%. Fi- 
nally, disc magnitudes have a standard deviation of 0.15 B 
mag. The largest residuals are mostly due to those objects 
that have different types between repeat observations (e.g. 
one observation results in a high B /T type 1 bugle-disc clas- 
sification, and the other observation results in a type 5 bulge 
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Figure 14. A summary of the steps involved in the logical filter for the Sersic+exponential catalogue. 



only classification in which cases only the bulge parameters 
are compared). 

Figure 1161 shows histograms of the distribution in dif- 
ferences between parameters from repeat observations using 
the de Vaucouleurs + exponential catalogue. Once again, the 
catalogues have been filtered according to the rules described 
in Section |5] These measurements also show an acceptable 
level of repeatability, although the scatter in the bulge differ- 
ences is larger than that found using the Sersic + exponen- 
tial catalogue. Bulge standard deviations are: AR e /R e — 
32.7%, At/e = 9.5%, and AM bu i ge = 0.29 mag. For discs, 
measurements are comparable to the Sersic + exponen- 
tial catalogue, with standard deviations of Ah/h = 6.7%, 
A cos(i) = 4.9% and A Af disc = 0.11 mag. 

Residuals for measured parameters for single- 
component Sersic fits are shown in Figure 1171 As would 
be expected this much simpler model produces much more 
repeatable results. In this case the standard deviations are 
A R e / R e = 3.7%, Af/t = 2.0%, AM = 0.04 mags, and 
A log(n) = 4.2%. 



7 FUNDAMENTAL BULGE AND DISC 
PROPERTIES 

To demonstrate the consistency and suitability of the 
separation into bulge and disc populations we recover 



two fundamental statistical results: the Kormendy relation 
jKormendvl Il977l) and the relation betw een disc central 
surface bright ness and disc scale- length dde Jond [l996b; 
lGrahamll200lh . In Figure [TH] the size-surface brightness re- 
lation is shown for all bulges with R e > 0.5T and Mb u i gc < 
— 17 mag, along with the relation of lKormendvl dl977Tl . The 
majority of red bulges follow the Kormendy relation well 
(especially for large values of B/T), however the smaller 
population of blue bulges do not. For reasons explained in 
iGraham fc Guzman! (12003!) , faint spheroids do not follow the 
Kormendy relation, having smaller R e , and fainter (^i e ). 

Figure 1191 shows the distribution in scale-length, ft, 
and central surface brightness, /io, for discs from the 
Sersic+exponential catalogue after filtering (see Section 
ED- Only galaxies with 1.678ft > 0.5r, and Mdi sc < -17 
mag (see Section are shown. The dashed line shows the 
AMbuige = —17 B mag limit. It is clear that larger discs have 
larger (i.e. fainter) values fo r their central surface brightne ss 
as shown by other authors ide Jondll996btlGrahanjl200lft . 

The right panel of Figure 1111 shows the distri- 
bution of bulges in the colour-log (n) plane from the 
Sersic+exponential catalogue (post-filtering). This distribu- 
tion c an be compared with that presented in iDriver et alJ 
( 2006) who consider total rather than component properties, 
and identify two distinct peaks. In Figure 1111 the bulges 
clearly form a single red, high-n peak, although there is a 
small shoulder of bluer, low-n objects. That the majority 
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Figure 15. Parameter comparison between repeat observations of galaxies. Bulge (top panel) and disc (lower panel) parameters are 
derived from Sersic-bulge/ exponential-disc filtered catalogues. The solid lines correspond to the mean difference between repeat mea- 
surements and the dashed lines show the standard deviation. In each panel the mean and standard deviation is noted. 



of bulges ide ntify with the red, high-n peak supports the 
conclusions of iDriver et alJ l)200(fl that the colour bimodal- 
ity reported in the galaxy population is best explained by 
two fundamental components: bulges and discs, rather than 
two different galaxy populations. Furthermore, this proba- 
bly reflects two dominant formation processes and epochs. 
The small number of blue, low-n bulges may also be evi- 
dence of a third population of pseudo-bulges resulting from 
secular evolution. The bimodal nature of the different com- 
ponents is again evident in Figure 1201 which shows the rela- 
tionship between the model fj,Q, and Mb for bulge, disc, and 
pseudo-bulge components. It is clear that there are two dis- 
tinct groupings with the high surface brightness bulges and 
the low surface brightness discs. The small number of blue, 
pseudo-bulges appear to have intermediate surface bright- 
nesses. 



8 CONCLUSIONS 

We have used GIM2D to perform 2D model fits (including 
bulge-disc decomposition) for 10095 galaxies with B M ac < 20 
mag from the Millennium Galaxy Catalogue. We initially 
produced three catalogues using: (1) a single- component 
Sersic model, (2) a Sersic+exponential bulge-disc model, 
and (3) an 7? 1 ' 4 +exponential bulge-disc model. We find that 
there is good agreement between GIM2D cent roids, mag- 
nitude s, and half-light radii, with those in the iLiske et alJ 
(2003) MGC catalogues. However, there are a significant 
number of objects that have incorrect mask images from 
SExtractor which had to be corrected manually. 

When a two-component model is used, we find that a 
significant fraction of galaxies have mathematically 'good' 
fits, which may not be the most appropriate or meaningful. 
The components may be inverted, cross each other once, 
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Figure 16. Comparison between repeat observations of galaxies for bulge (top panel) and disc (lower panel) parameters derived from 
R 1 ' 4 bulge + exponential disc filtered catalogues. The solid lines correspond to the mean difference between repeat measurements and 
the dashed lines show the standard deviation. In each panel the mean and standard deviation is noted. 



twice, or not at all. In addition, it is also important to stress 
that many galaxies are not two component systems, and sim- 
ply interpreting one component as 'bulge' and the other as 
'disc' would clearly be wrong. We identify eight distinct pro- 
file types based on which of the two components dominate 
at the galaxy centre, at fj, — 26 mag arcsec -2 , and whether 
the components cross. 

The majority (~ 50%) of objects are 'classical' type 1 
galaxies where the Sersic function dominates at the cen- 
tre (modelling a bulge), and the exponential function dom- 
inates the outer regions (modelling a disc). The remaining 
objects often require more careful interpretation. The differ- 
ent types are sometimes due to phenomena such as disc trun- 
cation, multiple nuclei, or strong irregularities, all of which 
lie beyond the scope of the simple GIM2D model. Using the 
'type' scheme we implement a strategy to ensure only gen- 
uine bulge-disc systems are fitted with a Sersic+exponential 



model, and replace all other objects with single-component 
Sersic fits. In most cases only single component systems 
end up with Sersic-only fits. A small number (< 6%) of ge- 
nuine bulge+disc systems are classified as type 3 or type 
4 (especially those with disc truncation), and will therefore 
be mis-classified as single-component systems. This 'logical 
filtering' process results in a fourth catalogue, where each 
galaxy is treated as either bulge+disc, bulge-only, or disc- 
only, and it is this catalogue that we suggest contains the 
most appropriate, physically meaningful modelling of the 
galaxy population. All five structural catalogues are made 
publically available with this paper (see Appendix A). 

To test the repeatability and accuracy of our measure- 
ments, we applied GIM2D twice to the objects that he in 
the overlap regions between the individual pointings that 
make up the MGC strip. We find that if size and absolute 
magnitude cuts are imposed, parameter measurements are 
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Figure 17. Comparison between repeat observations of galaxies for parameters derived from single component Sersic fits. The solid 
lines correspond to the mean difference between repeat measurements and the dashed lines show the standard deviation. In each panel 
the mean and standard deviation is noted. 




Figure 18. The (^ e ) - log(-Re) distribution for all MGC bulges. 
The red circles show objects with (u — r) corc 2.35, and the blue 
crosses denote objects with (u — ^)corc 

< 2.35. The dashed line 
sho ws the Mh„ i p- R = — 17 B mag limit, and the solid line shows 
the iKormendvl ll977l) relation. The red population follows the 
Kormendy relation well, with most of the scatter coming from 
objects with B/T < 0.5. The blue (pseudo-)bulges do not appear 
to fit the relation. 

repeatable at less than the 15% level. Although the cata- 
logues contain entries for all objects, we recommend impos- 
ing cuts in size (half-light radius > 0.5F) and component 
absolute magnitude (M < —17B mag) in order to obtain 
reliable measurements. To demonstrate the accuracy of the 
final, logically filtered catalogue we recover the fundamen- 



Figure 19. The fj,o — log(h) distribution for all reliable discs in 
the MGC (using the filtered Sersic+exponential catalogue). The 
dashed line shows the M = — 17 limit. Discs with larger scale- 
lengths have fainter central surface brightnesses. 



tal Kormendy and fj,o — log(Tl) relations for bulges and discs 
respectively. In a series of future papers, the catalogues pre- 
sented here will be used in detail to measure bulge and disc 
luminosity functions, bivariate brightness distributions, and 
size distributions. We are also using the catalogues to mea- 
sure the supermassive black hole mass function and to study 
the effects of disc opacity. 

Finally, we have used the catalogues to show that the 
galaxy colour-log(n) bimodality is due to the two-component 
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Figure 20. The distribution of galaxy component properties in 
the Mg — fio plane is shown by the dots. Red contours outline 
the distribution of bulge components, and blue contours outline 
the distribution of disc components. Two clear populations are 
evident. The magenta contours outline the blue ((u — r) coro < 
2.35) bulge components, which appear to lie in an intermediate 
position. The right panel shows histograms of the number of each 
type of component. The number of discs is boosted because they 
are more likely to meet the size criteria than bulges (or pseudo- 
bulges). 

nature of galaxies, rather than two distinct galaxy popu- 
lations. Luminous bulges generally occupy the red, high-n 
peak, and discs occupy the blue, low-n peak. There is also 
some evidence for a third population of blue, pseudo-bulges. 
The observed bimodality may be the end-result of two dis- 
tinct processes and timescales associated with bulge and disc 
formation. We conclude that routine bulge-disc decompo- 
sition is essential at all redshifts to fully understand the 
evolution of the luminous (Mb < — 17 mag) galaxy pop- 
ulation. Three major steps forward are now required: (1) 
high-resolution, deep near-IR survey data, (2) the expan- 
sion of bulge-disc decomposition software to accommodate 
real phenomena such as nuclei and disc truncation, and (3) 
deeper studies to characterise the dwarf (Mb > — 17 mag) 
population. 



ACKNOWLEDGEMENTS 

We thank Luc Simard for making his GIM2D code publi- 
cally available. The Millennium Galaxy Catalogue consists 
of imaging data from the Isaac Newton Telescope obtained 
through the ING Wide field Camera Survey Programme. 
The Isaac Newton Telescope is operated on the island of 
La Palma by the Isaac Newton Group in the Spanish Ob- 
servatorio del Roque de los Muchachos of the Instituto de 
Astrofisica de Canarias. Spectroscpic data comes from the 
Anglo- Australian Telescope, The ANU 2.3-m, the ESO New 
Technology Telescope, the Telescopio Nazionale Galileo, and 



the Gemini North Telescope. The survey has been supported 
through grants from the Particle Physics and Astronomy Re- 
search Council (UK), and the Australian Research Council 
(AUS). The data and data products are publically available 
from http://www.eso.org/~jliske/mgc/ or on request from 
J. Liske. 



REFERENCES 

Abazajian K., et al., 2003, AJ, 126, 2081 

Abraham R. G, van den Bergh S., Glazebrook K., Ellis 

R. S., Santiago B. X., Surma P., Griffiths R. E., 1996, 

ApJS, 107, 1 

Aguerri J. A. L., Iglesias-Paramo J., Vi'lchez J. M., Munoz- 

Tunon C, Sanchez- Janssen R., 2005, AJ, 130, 475 
Andredakis Y. C, Sanders R. H., 1994, MNRAS, 267, 283 
Baldry I. K., Glazebrook K., Brinkmann J., Ivezic Z., Lup- 
ton R. H., Nichol R. C, Szalay A. S., 2004, ApJ, 600, 
681 

Barden M., et al., 2005, ApJ, 635, 959 
Bertin E., Arnouts S., 1996, A&AS, 117, 393 
Blanton M. R., et al., 2003, ApJ, 594, 186 
Burbidge E. M., Burbidge G. R., Crampin D. J., 1964, ApJ, 
140, 1462 

Caon N., Capaccioli M., D'Onofrio M., 1993, MNRAS, 265, 
1013 

Capaccioli M., 1989, in World of Galaxies (Le Monde des 

Galaxies), pp. 208-227 
Christlein D., Zabludoff A. I., 2005, ApJ, 621, 201 
Ciotti L., Bertin G, 1999, A&A, 352, 447 
Cole S., Lacey C. G, Baugh C. M., Frenk C. S., 2000, 

MNRAS, 319, 168 
Colless M., et al., 2001, MNRAS, 328, 1039 
Conselice C. J., 2003, ApJS, 147, 1 

Cross N. J. G, Driver S. P., Liske J., Lemon D. J., Pea- 
cock J. A., Cole S., Norberg P., Sutherland W. J., 2004, 
MNRAS, 349, 576 

de Jong R. S., 1996a, A&AS, 118, 557 

— , 1996b, A&A, 313, 45 

de Jong R. S., Simard L., Davies R. L., Saglia R. P., 
Burstein D., Colless M., McMahan R., Wegner G, 2004, 
MNRAS, 355, 1155 

de Souza R. E., Gadotti D. A., dos Anjos S., 2004, ApJS, 
153, 411 

de Vaucouleurs G, 1948, Annales dAstrophysique, 11, 247 
— , 1959, Handbuch der Physik, 53, 275 
— , 1977, in Evolution of Galaxies and Stellar Populations, 
pp. 43-+ 

Doi M., Fukugita M., Okamura S., 1993, MNRAS, 264, 832 

Driver S. P., et al., 2006, MNRAS, 368, 414 

Driver S. P., Liske J., Cross N. J. G, De Propris R., Allen 

P. D., 2005, MNRAS, 360, 81 
Elmegreen B. G, Elmegreen D. M., Montenegro L., 1992, 

ApJS, 79, 37 
Fraser C. W., 1972, The Observatory, 92, 51 
Freeman K. C, 1970, ApJ, 160, 811 

Graham A., Lauer T. R., Colless M., Postman M., 1996, 

ApJ, 465, 534 
Graham A. W., 2001, MNRAS, 326, 543 
Graham A. W., de Blok W. J. G, 2001, ApJ, 556, 177 
Graham A. W., Driver S. P., 2005, PASA, 22, 118 



© 0000 RAS, MNRAS 000, 000-000 



20 Allen et al. 



Graham A. W., Guzman R., 2003, AJ, 125, 2936 
Graham A. W., Trujillo I., Caon N., 2001, AJ, 122, 1707 
Gutierrez C. M., Trujillo I., Aguerri J. A. L., Graham 

A. W., Caon N., 2004, ApJ, 602, 664 
Hubble E. P., 1926, ApJ, 64, 321 
Koo D. C, et al., 2005, ApJS, 157, 175 
Kormendy J., 1977, ApJ, 218, 333 
Lahav O., et al., 1995, Science, 267, 859 
Lahav O., Nairn A., Sodre L., Storrie-Lombardi M. C, 

1996, MNRAS, 283, 207 
Liske J., Lemon D. J., Driver S. P., Cross N. J. G., Couch 

W. J., 2003, MNRAS, 344, 307 
Mac Arthur L. A., Courtcau S., Holtzman J. A., 2003, ApJ, 

582, 689 

Madgwick D. S., et al., 2002, MNRAS, 333, 133 
Morgan W. W., 1958, PASP, 70, 364 
— , 1959, PASP, 71, 394 
— , 1962, ApJ, 135, 1 

Nairn A., et al., 1995a, MNRAS, 274, 1107 

Nairn A., Lahav O., Sodre L., Storrie-Lombardi M. C, 

1995b, MNRAS, 275, 567 
Odewahn S. C, Windhorst R. A., Driver S. P., Keel W. C, 

1996, ApJL, 472, L13+ 
Peng C. Y., Ho L. C, Impey C. D., Rix H., 2002, AJ, 124, 

266 

Poggianti B. M., 1997, A&AS, 122, 399 
Ravindranath S., et al., 2004, ApJL, 604, L9 
Saglia R. P., Bertschinger E., Baggley G., Burstein D., Col- 
less M., Davies R. L., McMahan R. K., Wegner G., 1997, 
ApJS, 109, 79 

Schade D., Lilly S. J., Crampton D., Hammer F., Le Fevre 

O., Tresse L., 1995, ApJL, 451, L1+ 
Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 

525 

Sersic J. L., 1968, Atlas de galaxias australes. Cordoba, 
Argentina: Observatorio Astronomico, 1968 

Simard L., Koo D. C, Faber S. M., Sarajedini V. L., Vogt 
N. P., Phillips A. C, Gebhardt K., Illingworth G. D., Wu 
K. L., 1999, ApJ, 519, 563 

Simard L., et al., 2002, ApJS, 142, 1 

Tamm A., Tenjes P., 2006, A&A, 449, 67 

Tasca L. A. M., White S. D. M., 2005, astro-ph/0507249 

Trujillo I., Aguerri J. A. L., 2004, MNRAS, 355, 82 

Trujillo I., et al., 2005, astro-ph/0504225 

van den Bergh S., 1976, ApJ, 206, 883 



© 0000 RAS, MNRAS 000, 000-000 



MGC: Bulge/Disc Decomposition of 10095 Galaxies 21 



APPENDIX A: CATALOGUES AND CATALOGUE PARAMETERS 

We provide five structural catalogues each containing 118 parameters for the 10 095 galaxies assembled from MGC-BRIGHT. 
The catalogues are available from the MGC website at http://www.eso.org/~jliske/mgc as: 



master. bsersic. cat 
master .dsersic.cat 
master. dve. cat 
master.se. cat 
master . logic . cat 



Catalogue containing Sersic-only profile fits with e < 0.7. 
Catalogue containing Sersic-only profile fits with i < 85.0. 
Catalogue containing de Vaucouleur plus exponential profile fits. 
Catalogue containing Sersic plus exponential profile fits. 

Bulge and disc catalogue after processing through the logical filter, containing Sersic+exponential bulge+disc 



The parameter file relevant for all five data files is: 



master. par - List of parameters included in the above five files. 
Full details for each column entry are given below. 



Column 1: (ID) MGC unique identification number as listed in NED numbers 00000-69999 represent original SExtractor 

detections and numbers 90000+ represent sources that have been rebuilt (see lLiske et al JI2003T) . 

Column 2: (B M gc) The SExtractor BEST magnitude correc ted for Galactic extinction via the ISchleeel et all (ll998T> dust 
maps, the photom etric system is defined in lLiske et alJ J2003) and conversions to various filters are listed in the appendix of 
ICross et al.U2004D . 

Column 3: ((/i e )) The apparent effective B M gc surface brightness inside the empirically measured seeing corrected half-light 
radius: (fi e ) = Smgc + 2.5 log 10 (27r[_R°] 2 ) (mag arcsec -2 ). 

Column 4 & 5: (RA & DEC) Right ascension and declination in J2000.0 (deg). 
Column 6: (CLASS) Classification parameter: l=galaxy, 8=star. 

Column 7: (R e ) The empirically measured half-light radius using the positional angle and ellipticity provided by SExtractor. 

The magnitudes are assumed total (arcsec). 

Column 8: (R°) Seeing corrected half-light radius ,i.e., R° = y/R% — 0.32r 2 fsee lDriver et al.ll2005h . (arcsec). 
Column 9: (Y) Full width at half maximum of the seei ng (arcsec ) . 

Column 10: (spec) Best spectral fitting template from iPoggiantil il997f) fsee lDriver et al]|2005t for fitting details). 
Column 11: (k z (B)) B M qc K-correction (mag). 
Column 12: (z) Best redshift for this galaxy. 

Column 13: (Q z ) Redshift quality flag: l=targeted but no redshift, 2= tentative redshift measurement, 3= reliable redshift, 
4 = definite redshift, 5= unequivocal redshift 9 = not targeted. 

Column 14: (msdss [Pet]) SDSS-DR1 extinction corrected Petrosian apparent magnitude (AB mag). 
Column 15: (<?sr>ss[Pet]) SDSS-DR1 extinction corrected Petrosian apparent magnitude (AB mag). 
Column 16: (rsusspet]) SDSS-DR1 extinction corrected Petrosian apparent magnitude (AB mag). 
Column 17: (isz>ss[Pet]) SDSS-DR1 extinction corrected Petrosian apparent magnitude (AB mag). 
Column 18: (zgosspet]) SDSS-DR1 extinction corrected Petrosian apparent magnitude (AB mag). 
Column 19: (ID) Duplicate of column 1. 

Column 20: (M USDSS [Pet]) SDSS-DR1 extinction corrected Petrosian Absolute 1 magnitude (AB mag). 
Column 21: (Mb moc [Kron]) MGC extinction corrected Kron Absolute 1 magnitude (B M gc mag). 
Column 22: (M gsDSS [Pet]) SDSS-DR1 extinction corrected Petrosian Absolute 1 magnitude (AB mag). 
Column 23: (M rsDSS [Pet]) SDSS-DR1 extinction corrected Petrosian Absolute 1 magnitude (AB mag). 
Column 24: (Mi SDSS [Pet]) SDSS-DR1 extinction corrected Petrosian Absolute 1 magnitude (AB mag). 
Column 25: (M ZSDSS [Pet]) SDSS-DR1 extinction corrected Petrosian Absolute 1 magnitude (AB mag). 

Column 26: (/j„ SDSS [Pet]) SDSS-DR1 extinction corrected Petrosian Absolute 1 effective surface brightness (AB mag arcsec -2 ). 
Column 27: (fiB M0C [Kron]) MGC extinction corrected Kron Absolute 1 effective surface brightness (B 

mgc rnag arcsec ) . 

Column 28: (fJ.g SDSS [Pet]) SDSS-DR1 extinction corrected Petrosian Absolute 1 effective surface brightness (AB mag arcsec" 2 ). 
Column 29: (fJ.r SDSS [Pet]) SDSS-DR1 extinction corrected Petrosian Absolute 1 effective surface brightness (AB mag arcsec" 2 ). 
Column 30: ()J-i SD ss P e t]) SDSS-DR1 extinction corrected Petrosian Absolute 1 effective surface brightness (AB mag arcsec -2 ). 
Column 31: (/iz SDSS [Pet]) SDSS-DR1 extinction corrected Petrosian Absolute 1 effective surface brightness (AB mag arcsec -2 ). 
Column 32: (M USDSS [PSF]) SDSS-DR1 extinction corrected PSF Absolute 1 magnitude (AB mag). 
Column 33: (M gsDSS [PSF]) SDSS-DR1 extinction corrected PSF Absolute 1 magnitude (AB mag). 
Column 34: (M rsDSS [PSF]) SDSS-DR1 extinction corrected PSF Absolute 1 magnitude (AB mag). 
Column 35: (Mi SDSS [PSF]) SDSS-DR1 extinction corrected PSF Absolute 1 magnitude (AB mag). 
Column 36: (M ZSDSS [PSF]) SDSS-DR1 extinction corrected PSF Absolute 1 magnitude (AB mag). 

Column 37: (8i) Angular size of lh -1 pc at object redshift: 8i = 3600 tan -1 [-jggg^-], a value of 0.0 implies no redshift (arcsec). 
Column 38: (d p ) Proper (co- moving) distance to the object (ft -1 Mpc). 

Column 39: (c type) Continuum type: 1=E1 15 Gyr, 2= Sa 7.4 Gyr, 3= Sc 2.2 Gyr fsee lDriver et alJl2006^ . 

Column 40: (m type) Eyeball morphological type: = not classified, 1 = E/S0, 2— Sabc, 3=Sd/Irr (all Bmgo < 19 mag have 

been classified). 
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ID) Duplicate of column 1. 

2dFGRS No.) Matched 2dFGRS serial n umber for this object (000000 if no match). 

•q) 2dFGRS r\ parameter if matched (see iMadewick et al.l l2002:. -99.9 if no match). 

ID) Duplicate of column 1. 

L_TOT) Total flux (digital units). 

L.TOT-) Total flux error (-). 

L.TOT+) Total flux error (+). 

BULGE_FRAC) Bulge fraction (0=pure disc). 

BULGE_FRAC-) Bulge fraction error (-). 

BULGE_FRAC+) Bulge fraction error (+). 

BULGE_RE) Bulge effective radius (pixels). 

BULGE_RE-) Bulge effective radius error (-). 

BULGE.RE+) Bulge effective radius error (+). 

BULGE_E) Bulge ellipticity. 

BULGE_E-) Bulge ellipticity error (-). 

BULGE_E+) Bulge ellipticity error (+). 

BULGE_PA) Bulge position angle. 

BULGE_PA-) Bulge position angle error (-). 

BULGE_PA+) Bulge position angle error (+). 

R_D) Exponential disc scale length (pixels). 

R_D-) Exponential disc scale length error (-). 

R_D+) Exponential disc scale length error (+). 

DISC J) Disc inclination (0=face-on). 

DISCJ-) Disc inclination error (-). 

DISCJ+) Disc inclination error (+). 

DISC_PA) Disc position angle. 

DISC_PA-) Disc position angle error (-). 

DISC_PA+) Disc position angle error (+). 

X_OFF) X offset of galaxy centre (pixels). 

X_OFF) X offset error (-). 

X_OFF) X offset error (+). 

Y_OFF) Y offset of galaxy centre (pixels). 

Y_OFF) Y offset error (-). 

Y_OFF) Y offset error (+). 

BACK) Background level (digital units). 

BACK-) Background level error (-). 

BACK+) Background level error (+). 

n) Sersic index (=4 for de Vaucouleurs profile). 

n-) Sersic index error (-). 

n+) Sersic index error (+). 

CHI) x 2 of GIM2D fit. 

RHALF) GIM2D HLR (pixels). 

CI) Concentration index alpha=l. 

Al) Asymmetry index (background corrected) alpha=l. 
Bl) Background correction applied to Al. 
C2) Concentration index alpha=2. 

A2) Asymmetry index (background corrected) alpha=2. 
B2) Background correction applied to A2. 
C3) Concentration index alpha=3. 

A3) Asymmetry index (background corrected) alpha=3. 
B3) Background correction applied to A3. 
C4) Concentration index alpha=3. 

A4) Asymmetry index (background corrected) alpha=3. 
B4) Background correction applied to A3. 

Mb moc [GIM2D]) Total extinction corrected absolute 2 magnitude derived from L_TOT (mag). 

M Bmgc (Bulge) [GIM2D]) Total extinction corrected absolute 2 magnitude derived from BULGE_FRAC*L_TOT 

M Bmgc (Disc) [GIM2D]) Total extinction corrected absolute 2 magnitude derived from (l-BULGE_FRAC)*L_TOT 



R e (Bulge) [GIM2D] ) Half-light radius of bulge component: R e = O.333BULGE_E/0 (h^ 1 
Qo(Disc)[GIM2D]) Scale-length of disc component: a D = 0.333R_D/6» (ft -1 kpc). 



kpc). 
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Column 100: (no (Bulge)) Absolute central surface brightness: fi (Bulge) = fi e (Bulge) — 1.08576 fsee lGraham fc Driveif 2005. 
eqn. 7) (mag arcsec -2 ). 

Col umn 101: (u e (Bulge)) A bsolute effective surface brightness at R e (Bulge): /^(Bulge) = (/x e ) (Bulge) + 2.5 log 10 [^ r r(2n)] 
fsee lGraham fc Driv er 2005, cqn. 9) (mag arcsec -2 ). 

Column 102: ((/^ e ) (Bulge) ) Absolute effective sur face brightness of bulge component: (/x e ) (Bulge) = Ms MG0 (Bulge) + 
2.51og 10 (27nR 2 ) +36.57 fsee lGraham fc Drive JI200.4 eqn. 12) (mag arcsec" 2 ). 

Column 103: (fio (Disc)) Absolute central surface brightness: /i (Disc) = /i e (Disc) — 1.08576 (see iGraham fc Driverl l2005. 
eqn. 7) (mag arcsec -2 ). 

Column 104: (lie. (Disc) ) Absolute effective surface brightness at 1.678an: /i e (Disc) = (/i c )(Disc) + 2.5 log 10 [p^-r(2n)] (see 
IGraham fc Driverll2005l eqn. 9) (mag arcsec -2 ). 

Column 105: ((/i e ) (Disc)) Absol ute effective surface b rightness of disc component: (/i e )(Disc) = Mb mg0 (Disc) + 
2.51og 10 (27r(1.678a D ) 2 ) +36.57 f see IGraham fc Driverlbool eqn. 12) (mag arcsec -2 ). 

Column 106: (R x ) Radius at which the surface brightness of the disc equals the surface brightness of the bulge, i.e., 
< M e)(Bulge) + 1.08576( fle( g- lge) )^ = (/ie)(Disc) + 1.0857(|j) fseelGraham fc Driverl!^ eqn. 17) (arcsec). 

Column 107: (Hx) Surface brightness at R x (mag arcsec— 2). 
Column 108: (p type) Profile type as defined in Section [5] 

Column 109: (b) the value which satisfies T(2n) = 27(2n, b) fsee lGraham fc Driver! 1 2005). 
Column 110: (ID) Duplicate of column 1. 

Column 111: ((u — r) g ) Rest frame global colour: M u (Pet) — M r (Pet) (AB mag). 
Column 112: ((u - r) c ) Rest frame core/bulge colour: M„(PSF) - M r (PSF) (AB mag). 

10 -0.4(u-r) a _ B 1() -0.4(u-r) c 

Column 113: ((u — r)d) Rest frame outer/disc colour: (u — r)d = — 2.51og 10 [ % ] (AB mag). 

1— -y 

Column 114: (ID) Duplicate of column 1. 

Column 115: (FLAGS) SExtractor flags from original detection. 
Column 116: (ID) Duplicate of column 1. 

Column 117: (S) Galaxy significan ce: i.e., £%9^ ( seelDriver et alJl200ll2006l) . 
Column 118: (AS) Error in S (see iDriver et al.ll2005l Eoolj) . 

Notes: 

1 Absolute magnitudes are derived using Qm = 0.3, Q.a = 0.7, Ho = 100 km s -1 Mpc -1 , the appropriate k-correction derived 
for each filter and an evolutionary correction of the form I/ z= o.o = L z (l + z)~^ where /3 = 1.5,0.75,0.5,0.375,0.3,0.2 for 
u, B MGC , g,r,i, z respectively. Galaxies without redshifts have values of -99.9 in these columns. 

2 GIM2D absolute magnitudes are derived using: M T g D = Z p - 2.5 log 10 (X T g D L_TOT) - (B - Mb) - A b where X = 1, 
BULGE_FRAC or (1-BULGE_FRAC) for Total (T), Bulge (B) or Disc (D) respectively (null values are set to 99.99). 
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